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CYCLOSTATIONARY IN THE TIME VARIABLE UNIVERSE 
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ABSTRACT 

Cyclostationary processes are those signals whose have vary almost periodically in statistics. It can 
give rise to random data whose statistical characteristics vary periodically with time although these 
processes not periodic functions of time. Intermittent pulsar is a special type in pulsar astronomy 
which have period but not a continuum. The Rotating RAdio Transients (RRATs) represent a 
previously unknown population of bursting neutron stars. Cyclical period changes of variables star 
also can be thought as cyclostationary which are several classes of close binary systems. Quasi- Periodic 
Oscillations (QPOs) refer to the way the X-ray light from an astronomical object flickers about certain 
frequencies in high-energy (X-ray) astronomy. I think that all above phenomenon is cyclostationary 
process. I describe the signal processing of cyclostationary, then discussed that the relation between 
it and intermittent pulsar, RRATs, cyclical period changes of variables star and QPOs, and give the 
perspective of finding the cyclostationary source in the transient universe. 

Subject headings: Cyclostationary — methods: data analysis: Intermittent pulsar, RRATs, Cyclical 
period changes of variables star, QPOs 



1. INTRODUCTION 

It was first mentioned by the design of synchroniza- 
tion algorithms for communications syst ems that cy- 
clostationary process ()Serpedin et al.l [2005h . Many pro- 
cesses encountered in nature arise from periodic phe- 
nomena. For example, in telecommunications, telemetry, 
radar, and sonar applications, periodicity is due to mod- 
ulation, sampling, multiplexing, and coding operations. 
In mechanics it is due, for example, to gear rotation. 
In econometrics, it is due to seasonality; and in atmo- 
spheri c science it is due to rotation and revolution of the 
earth (|Gardner et al.l200llSerpedin et alj2005HGardnerl 
I1994D . Cyclostationary processes are named in multiple 
different ways such as periodically correlated, periodi- 
cally nonstationary, periodically nonstationary or cyclic 
correlated processes in literature (jHurdl 1 1989L 119971 ). 

In astronomy, now just in removing cyclostationary ra- 
dio frequ ency interferences (RFI) from radio astronom - 
ical data (|Leshem et al.l 120001: iBretteil fc Weber! 120051) 
and the galactic w hite-dwarfs backgro und that will be 
observed by LISA(Ed lund et al.l 120051 ) cyclostationary 
process be mentioned. Some radio pulsar have been 
discover ed the phenomenon of intermittent period in 
recently (|Kramer et al.l I2006T) . The nulling of the pul- 
sar be thought a characteristic of the older pulsars. In 
X-ray pulsar, the phenomenon of intermittent far-flung 
exist for the effect of companion. It also show almost 
period in time series. The Rotating RAdio Transients 
(RRATs) demonstrate a pre viously unknown popula ,- 
tion of bursting neutron stars(|McLaughl m et al.l 12006). 
Cyclical period changes of variables star that includ- 
ing Algol, WUrsaeMajoris, and RS Canum Venaticorum 
systems and the cataclysmic variables, and RR Lyrae 
star et al. also can be thought as cyclostationary. The 
fastest variability components in X-ray binaries are the 
kilohertz quasi-periodic oscillations (kHz QPOs), which 
occur in a wide variety of low magnetic-field neutron 
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star svstems ljvan der Klisl I2005T ). The signal process of 
the above phenomenon, in a statistical sense, as a peri- 
odic function of time. These kind of processes have been 
studied for many years, and are usu ally referred to as cy- 
clostationary random processes (see ([Gardner et al.ll2006t 
iSerpedin et alJl2005l !Gardnerlll994l ) for a comprehensive 
overview of the subject and for more references). 

In what follows we will briefly summarize the prop- 
erties of cyclostationary processes, then discuss the cy- 
clostationary signal process in some astronomical phe- 
nomenon, in the end give the perspective of finding cy- 
clostationary processes in the time variable universe. 



2. THE BASIC OF CYCLOSTATIONARY PROCESSES 

The continuous stochastic process X(t) having finite 
second order moments is said to be cyclostationary with 
period T if the following expectation values 



E[X(t))=m{t) = m(t + T), 
E[X(t')X{t)\ = C(t', t) = C(t' +T,t- 



T) 



(1) 
(2) 



are periodic functions of period T, for every (t',t) 6 
R x R. We will assume m(t) = for simplicity 
now ljOardner fc Franksl[l975l: IHiir dlll989L fl997h . The im- 
portant special case of cyclostationary signals are those 
which exhibit cyclostationary in second-order statistics 
(e.g., the autocorrelation function). It is called wide- 
sense cyclostationary signals, and is analogous to wide- 
sense stationary processes. The exact definition differs 
depending on whether the signal is treate d as a stochast ic 
process or as a deterministic time series (|Gardnerlll978l ). 

If X{t) is cyclostationary, then the function B(t, r) = 
C(t + r, t) for a given r e R is periodic with period T, 
and it can be represented by the following Fourier series 



B(t,r) 



T — — OO 



B r {r)e 



(3) 



2 



where the functions B r (r) are given by 

B r (r) = ^ £ B{t,r)e- i2 ^ dt . (4) 

The Fourier transforms g r (f) of B r (r) are the so called 
"cycl ic spectra" of the cyclostationary process X(t) ( Hurd 
[19891) 

/OO 
B r { T )e~ a ^ T dr . (5) 
-oo 

If a cyclostationary process is real, the following relation- 
ships between the cyclic spectra hold 

B^{t)=B*{t) , (6) 
9-r(-f)=9*Af) , (7) 

where the symbol * means complex conjugation. This 
implies that, for a real cyclostationary process, the cyclic 
spectra with r > contain all the information needed to 
characterize the process itself. 

The function cr 2 (r) = B(Q,t) is the variance of the 
cyclostationary process X(t), and it can be written as a 
Fourier decomposition as a consequence of Eq. ((4]) 

oo 

a 2 (r)= J2 H r e^9, (8) 

r— — oo 

where H r = B r (0) are harmonics of the variance a 2 . 
From Eq. © it follows that H_ r = H*. 

For a discrete, finite, real time series X t , t = 1, . . . , M 
the cyclic spectra can be estimated by generalizing stan- 
dard methods of spectrum estimation used with station- 
ary processes. Assuming again the mean value of the 
time series X t to be zero, the cyclic autocorrelation se- 
quences are defined as 

N-\l\ 

s i =J/J2 X t X t+w e- ,3sl ^ . (9) 

The cyclic autocorrelations are asymptotically (i.e. 
for N — > oo) u nbiased estimators of the functions 
Br-ir) ( Hurdl 1 19891 ) . The Fourier transforms of the cyclic 
autocorrelation sequences s\ are estimators of the cyclic 
spectra g r {f)- These estimators are asymptotically un- 
biased, and are called "inconsistent estimators" of the 
cyclic spectra, i.e. their variances do not tend t o zero 
asym ptotically. In the case of Gaussian processes (|Hurdl 
1989) consistent estimators can be obtained by first ap- 
plying a lag window to the cyclic autocorrelation and 
then perform a Fourier transform. 

The alternative procedure for identifying consistent es- 
timators of the cyclic spectra is to first take the Fourier 
transform (|Gardner et al.ll2006l : IGardnerlll994l) . X(f), of 
the time series X (t) 

AT 

X(f) = Xte-^^ (10) 
and then estimate the cyclic periodograms g r (f) 

9r(f) = U> f T — • (11) 



By finally smoothing the cyclic periodograms, consistent 
estimators of the spectra g r (f) are then obtained. The 
estimators of the harmonics H r of the variance a 2 of a 
cyclostationary random process can be obtained by first 
forming a sample variance of the time series X t . The 
sample variance is obtained by dividing the time series 
X t into contiguous segments of length tq such that tq 
is much smaller than the period T of the cyclostationary 
process, and by calculating the variance erf over each seg- 
ment. Estimators of the harmonics are obtained either 
by Fourier analyzing the series a 2 or by making a least 
square fit to a 2 with the appropriate number of harmon- 
ics. Note that the definitions of (i) zero order (r = 0) 
cyclic autocorrelation, (ii) periodogram, and (iii) zero or- 
der harmonic of the variance, coincide with those usually 
adopted for stationary random processes. Thus, even 
though a cyclostationary time series is not stationary, 
the ordinary spectral analysis can be used for obtaining 
the zero order spectra. Note, however, that cyclostation- 
ary random processes provide more spectral information 
about the time series they are associate d with due to 
the e xistence of cyclic spectra with r > O ljEdlund et al.1 
l2005h . 

For stationary and cyclostationary time-series, there 
is two alternative philosophical frameworks for the two 
problems of estimating the time- invariant or time- variant 
autocorrelation function and its Fourier transform. One 
is based on the stochastic process model, and the other is 
based on the nonstochastic time-series model. Gardner, 
W.A. compared it, and explained that results on estima- 
tor bias and variance for these two problems couched 
within the stochastic process framework have analogs 
within the nonstochastic framework. The bias and vari- 
ance results for cyclostationary ti me-series that a re avail- 
able within these two frameworks (| Gardner! 1 1 9 9 ll ) . 

As an important and practical application, assuming 
consider a time series yt consisting of the sum of a sta- 
tionary random process, nt, and a cyclostationary one 
X t (i.e. yt — n t + X t ). Let the variance of the stationary 
time series n t be v 2 and its spectral density be £(/). It 
is easy to see that the resulting process is also cyclosta- 
tionary. If the two processes are uncorrelated, then the 
zero order harmonic of the variance of the combined 
processes is equal to 

Z 2 = v 2 + a 2 , (12) 
and the zero order spectrum, Go(f), of yt is 

Go(f) = £(/)+ go (/) ■ (13) 

The harmonics of the variance as well as the cyclic 
spectra of yt with r > coincide instead with those 
of X t . In other words, the harmonics of the variance 
and the cyclic spectra of the process y t with r > 
contain information only about the cyclostationary pro- 
cess X t , and are not "contam inated" by the stationary 
processn t (jEdlund et al.l l2005| ) . 

3. CYCLOSTATIONARY IN PULSAR, QPOS AND 
VARIABLES STAR 

Recently, Kramer, M. et al. discovered one class of 
neutron stars that are seemingly ordinary radio pulsars, 
but which are only active for some short time and in a 
quasi-periodic fashion. So they call these "Intermittent 



3 



Pulsars" . These pulsar that is only periodically active. 
It appears as a normal pulsar for about a week and then 
" switches off' for about one month before emitting pulses 
again. The pulsar, called PSR B1931+24, is unique in 
this behaviour and affords astronomers an opportunity to 
compare its quiet and active phases. Most surprisingly, 
the pulsar rotation sl ows down faster wh en the pulsar is 
on than when it is off (|Kramer etali r2006). In the Fi gure 
[U part a) show a typical sequence of observations cover- 
ing a 20-month interval is indicated by the black lines. It 
shows respectively the times of observation and the times 
when PSR B1931+24 was on. It is clear that the pul- 
sar is not visible for ~80% of the time, part b) give the 
appearance of the pulsar is quasi-periodic nature, demon- 
strated by the power spectrum of the intensity obtained 
from the Fourier Transform of the autocorrelation func- 
tion of the mean pulse flux density obtained over the 
same 20-month interval, part c) is histograms of the du- 
rations of the on (solid) and off (hatched) phases. In 
off phases, integration over several weeks shows that any 
pulsed signal has a mean flux den sity of less than 2 /it Jy 
at 1400 MHz rtKramer et al.ll20M ).. 

McLaughlin, M. A. et al. discovered a class of Rotating 
Radio Transients(RRATS) which are identified as rotat- 
ing neutron stars that send out very short flashes of ra- 
dio light. These flashes are very short and very rare: one 
hundredth of a second long, the total time the objects 
are visible amounts to only about one tenth of a second 
per day. The isolated flashes last for between 2 and 30 
milliseconds. In between, for times ranging from 4 min- 
utes to 3 hours, the new stars are silent. They current 
estimates suggest that these objects are four times more 
common in the Galaxy than radio pulsars. In the Fig- 
ure [21 from top to bottom, it is show that the original 
detections of J1317-5759, J144 3-60 and J1826-14 in th e 
Parkes Multibeam Survey data ljMcLaughlin et al.l l2006) . 
It is show that RRATs still have period in short time 
scale. They have identified peri odicities in the range of 
0.4 to 7 s in 10 of the II obiects (|Lvnell2007f) . 

Camilo, F. et al. show that XTE J18102197 emits 
bright, narrow, highly linearly polarized radio pulses, 
observed at every rotation, thereby establishing that 
magnetars can be radio pulsars. There is no evidence 
of radio emission before the 2003 X-ray outburst (un- 
like ordinary pulsars, which emit radio pulses all the 
time), and the flux varies from day to day. The flux 
at all r adio frequencies is a pproximately equal-and at > 
20GHz l|Camilo et al.l[20M) . 

In wide-sense, if one pulsar have nulling or giant pulse 
phenomenon, it will not have strict period function with 
time, just in statistics have periodic. Another pulsar cy- 
clostationary process is Shabanova, T.V. find that the 
nature of the observed cyclical changes in the timing 
residuals from PSR B1642 03 is a continuous genera- 
tion of peculiar glitches whose amplitudes are modulated 
by a periodic large-scale sawtooth-like function. As the 
modulation function is periodical, the picture of cyclical 
timing residuals will be exactl y repeated in eac h modu- 
lation period or every 60 years (jShabanovall2009f ). 

QPOs were first id e ntified in white dwarf 
svstems( |van der Klis et all I1985D a nd then in neu- 
tron star systems (Jvan der Klisl 12005). Two QPO peaks 
(the 'twin peaks') occur in the power spectrum of the 
X-ray flux variations. They move up and down in fre- 



quency together in the 300- 1200Hz range in correlation 
with source state and often, luminosity. The typically 
300-Hz peak separation usually decreases by a few 
ten s of Hz when bot h peaks move up by hundreds of 
Hz (|van der Klisl 120051 ) . In the Figure 02 show that Keck 
II spectroscopy of optical mHz quasi-periodic oscillations 
(QPOs) in t he light curve of th e X-ray pulsar binary 
Hercules X-l (|Q'Brien et al.|[200l . 

Year- to decade-long cyclic orbital period changes have 
been observed in several classes of close binary systems, 
including Algol, W Ursae Majoris, and RS Canum Ve- 
naticorum systems and the cataclysmic variables. The 
origin of these changes is unknown, but mass loss, apsi- 
dal motion, magnetic activity, and the presence of a third 
body have all been proposed. In the Figured show that 
cyclical period changes in the dwarf novae V2051 Oph. 

The other two statistics periodic phenomenon also 
relate with stellar. One is Hallinan, G. et al. de- 
tected periodic bursts of extremely bright, circularly 
polari zed, coherent radio emission from the ultracool 
dwarf (Ha llinan et al.ll2007l ). Another is Double Periodic 
Variables (DPVs) that are blue stars characterized by a 
short periodicity (1-16 days) and a long periodicity (50- 
600 days) in their light curves. They were discovered in 
the Magellanic Clouds after a search for Be stars in the 
OGLE variable star catalog! Mennic kent et al. l l2003f ) 

4. DISCUSSION AND CONCLUSIONS 

Signal detection techniques designed for cyclostation- 
ary signals take account of the periodicity or almost pe- 
riodicity of the signal autocorrelation function. Single- 
cycle and multicycle detecto rs exploit one or mult iple cy- 
cle frequencies, respectivelv (|Gardner et al.ll2006[ ). Some 
search- efficient methods of det ection of cyclostationary 
signals (|Yeung fe Gardner! 1 19961 ) an d higher-order cyclo- 
statio nary for weak-signal detection ()Spooner fe Gardner! 
Il992f ) be developed. These will benefit of search tran- 
sient source in the time variable universe, especially that 
have periodic in statistics. 

Although no periodicities were detected in any of the 
sources using standard Fourier or folding methods, for 
ten of the sources (RRATs) McLaughlin, M.A. et al. iden- 
tify a periodicity from the arrival times of the individual 
bursts that used search technique s similar to those de- 
scribed in (jCordes fe M cLaughlin! I2003T ). In short, the 
35-minute time series were dedispersed for a number of 
trial values of DM. The time series were smoothed by 
convolution with boxcars of various widths to increase 
sensitivity to broadened pulses, with a maximum box- 
car width of 32 ms. Because the optimal sensitivity is 
achieved when the smoothing window width equals the 
burst width, our sensitivity is lower for burst durations 
greater than 32 ms. Each of these time series was then 
searched for any bursts above a threshold of five standard 
deviations, computed by calculating a running mean and 
root-mean-square deviation of the noisy time series. All 
bursts detected above a 5-cr threshold are plotted as cir- 
cles, with size proportional to the signal-to-noise ratio of 
the detected burst. The abcissa shows arrival time while 
the ordinate shows the DM. Because of their finite width, 
intense bursts are detected at multiple DMs and result 
in vertical broadening of the features. Bursts which are 
strongest at zero DM and therefore likely to be impul- 
sive terrestrial interference are not shown. In general 
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these were easily identified by t heir detection in multipl e 
beams of the 13-beam receiver (|McLaughlin et al.ll2006f) . 
More recently, Deneva, J.S. et al. use the above method 
and a friends-of-friends algorithm perform the ongoing 
Arecibo Pulsar ALFA (PALFA) survey of the Gala ctic 
planethen discover seven objects (|Deneva 

The discovery of RRATs increases the current Galactic 
population estimates of radio pulsar by at least several 
times. It seems also that there will be many candidates 
for whom it will be impractical or impossible to follow 
up at present with current observing facilities. These will 
require followup with instruments like LOFAR, FAST or 
the SKA. Keane, E.F. et al. note that these instruments 
will produce extraordinarily large volumes of data so that 



searching for transient RRAT-like sources will necessi- 
tate the development of automa ted algorithms whi ch will 
use the steps as outlined above ljKeane et af1 l2"009). So I 
think we should use the detection methods of cyclosta- 
tionary process integrate with DM search, that will be 
one automated and effective algorithms to seek RRATs. 

If we do similar things in find the other astronomi- 
cal cyclostationary sources that include Intermittent Pul- 
sars, QPOs et al., it will lead to discover more unusual 
astronomical phenomenon. 



DJ thanks 
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Fig. 1. — Time variation of the radio emission of PSR B1931+24. During the on phases, the pulsar is easy to detect and has the stable 
long-ter m intrinsic flux den sity associated with most normal pulsars. Since 1998, the pulsar has been observed as frequently as twice a 
dav.Csee lKramer et alj|20()6ft 
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Fig. 3. — Left panel, lightcurv e of the optical varia bility in Hercules X-l for the entire data-set. Right panel, lightcurve of a subset of the 
data clearly showing the QPO. (O'Brien et al. 2001) 
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Fig. 4. — The dashed and dot-dashed lines in the lower panel sho w the best- fit 11 yr c ycle period sinusoidal ephemeris, respectively, for 
the data in the first and the second halves of the time interval, (see Baptista ct al. 2003) 



